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CAA For J et Noise Physics - Reda Mankbadi, Embry-Riddle Aeronautical 
University as summarized by Eugene Krejsa. 

Dr. Mankbadi summarized recent CAA results. Examples of the effect of various 
boundary condition schemes on the computed acoustic field, for a point source in a 
uniform flow, were shown. Solutions showing the impact of inflow excitations on the 
result were also shown. Results from a large eddy simulation, using a fourth-order 
MacCormack scheme with a Smagorinsky sub-grid turbulence model, were shown 
for a Mach 2.1 unheated jet. The results showed that the results were free from 
spurious modes. Results were shown for a Mach 1.4 jet using LES in the near field 
and the Kirchhoff method for the far field. Predicted flow field characteristics were 
shown to be in good agreement with data and predicted far field directivities were 
shown to be in qualitative agree with experimental measurements. 

Dr. Mankbadi also presented results using linearized Euler equations. Comparison 
of predicted directivities agreed well with measurements for supersonic jet. 
Agreement was not as good for jet with Mach numbers less than 0.9. Results from 
very large eddy simulation were also presented. Dr. Mankbadi concluded his 
presentation with observations that CAA can provide: 

• Prediction of sound propagation. 

• Physics of the very large flow structures. 

• Correlation for other approaches, such as MGB. 

• Numerical experiments with various ideas for identification and 
control of the noise sources. 

There are still many challenges: 

• LES is CPU intensive 

• Shock-acoustic wave interaction is difficult 

• Coupling to the engine requires curvilinear meshes, grid generations, 
and a much faster code. 
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OUTLINE 

1. Numerical Issues 

2. LES/VLES for Jet Noise 

3. Physics-based Methods for Jet Noise Predictions 
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Numerical Issues 


• Large Computational Domain (60 D x 40 D) 

• Three-Dimensional Problem 

• Large Range of Length Scales 

• Acoustics: Acoustic Wavelength 

• Large Scale Structures: Jet Diameter 

• Shear-Layer Scales: Momentum Thickness 

• Mean Flow 

• High-accuracy numerical scheme is needed to resolve these scales in 

large domain 

• Dispersion and dissipation errors must be small for amplitude and 

phase accuracy 

• For linear problems, several schemes are available. 

• For mixing noise (nonlinear), few schemes are suitable: DRP, High- 

accuracy MacConnack, Compact, Prefactored Compact 

• Boundary Conditions are still the major difficulty 
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Technical 

Evaluation of Boundary Conditions on Benchmark Cascade Case 

• In this benchmark problem, a periodic vortical gust impinges on a flat plate cascade. The dose proximity of the 

computational boundaries makes this problem extremely difficult The results of several boundary conditions are shown. 
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EFFECT OF BOUNDARY TREATMENT ON ACCOUSTIC RADIATION 
POINT SOURCE IN A UNIFORM FLOW 
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Exact Solution 


Giles Boundary Condition 




Tam and Webb Boundary Condition 


Thompson Boundary Condition 
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AXISYMETRIC COMPUTATION FOR ACOUSTIC 
§ RADIATION FROM AXI SYMMETRIC LARGE 
SCALE STRUCTURE 
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Boundary Conditions 

Acoustic Radiation 
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Acoustic Radiation 


Mj = 2.1 
u^Uj = 0 

Re = 70000 
T, = 294°k 
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Instantaneous Pressure Distribution For M=2. 1 Axisymmetric Jet 


Single Frequency, e=0.04 


Single Frequency, e=0.001 



Bi- Modal Excitation, e=0.04 


Random Frequency, e=0.04 
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INFLOW EXCTTAITONS 


• Single Frequency: 

[u 1 v' p’ p] = e5R (r) e* ~ Wt) } 
e = input excitation level, a = eigenvalue. 

O(r) is the corresponding eigenfunctions $ (r) = [(i(r) $(r) p(r) p(r 


• Bi-Modal Excitation: 


T 1 at r A / i (ax -cot) { *f a s x ®s t + 0l 

[u' v’ p‘ pj = e9t{0(r)e } + e9t1 <E> s (r)e ^ 'f 


where the subscript s represents subharmonic frequency, 
p = initial phase difference between the fundamental and the subharmonic. 


Random Disturbance: 

[u' v- p' p] = eA(t)exp[-ln(2)(LzhMj 2 


A(t) is a random function in time, created by random function generator. 

b(0) and h(0) are the half-width and potential core radius of the shear layer at jet inflow. 


AXISYMMETRIC SIMULATIONS 


A fourth-order scheme is appropriate 

Clean solution can be achieved via careful attention to boundary 
treatment 

The solution is dependent on the inflow disturbances 
The wavelike nature of the large-scale structure 



LARGE- EDDDY SIMULATIONS 


Fourth-Order MacCormack 
Smagronisky Model 
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AXISYMMETRIC JET, unheated, fully expanded, 
M=2.1 . St=0.2 & 0.4 

COMPUTATIONAL DETAILS: 

Uniform grid spacing in the axial & Azimuthal 
directions. 

Stretched grid in the radial direction 
The domain: 5<x/R<70, 0<r/R< 32 
theta O to pi 
CFL=0.5 

25 points per wavelength 

675 time steps per cycle, 100 time steps per 
cycle 

AXIAL (196) radial (385), x=391 r=300 

x=1 96, r=385, n=12 

RESULTS: 

Though viscous and damping effects are 
damping the solution is free from spurious 
modes. 
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M=2.1 Supersonic Jet Excited by First Helical Mode 



Axial Velocity 



Vorticity 
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Technical Hig hli ght 

* f 

M=1.4 supersonic round jet excited by multiple frequencies at nozzle exit: 


Large-Eddy Simulation in ; the near field 



Mach Number 



Kirchhoff Method in the far field. 





Enhanced mixing and the breakdown of 
organized structures downstream of jet 
are observed. The far field noise directivity 
shows qualitative agreement with the 
experimental measurement 
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3D LES 

• Direct, 3D predictions is expensive 

• Need to split the problem into a source regime in the near field and a 
propagation regime and use physically based approximation * 
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CAN LEE APPROXIMATE BOTH THE NEAR AND FAR FIELD ? 
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Snapshot of Oscillating Pressure Field 
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Frequency Spectra in the Shear Layer 
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Frequency Spectra in the Far Field 



(x/D =10; r/D= 11.75) 
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Far Field Radial Decay 



R/D 
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Noise levels for Axisymmetric Mode 


(Symbols: current numerical calculation. Lines: Analytical solution) 
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Overall Comparison (St = 0.2) 
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Directivity of Jet Noise (R = 24) 
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Directivity of Jet Noise (R = 24) 
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SOUND PROPGATION- LEE 


Assuming the mean flow is assumed to be given, Linear Euler Equations 

(LEE) can be used for prediction of sound propagation and for 

approximate prediction of the sound source 

PHYSICS 

• Errors may amplify more in the linear simulations! 

• Need to make sure that no spurious modes are generated 

• The solution should decay as 1/r in the far field 

• Excellent agreement with Tam & Burton asymptotic analysis for 
instability waves 

• Reasonable agreement with Trout & McLaughlin M=2.1 controlled 
experiment 

• LEE can provide reasonable accurate approximate prediction for 
the large scale structure as low as M=0,9 
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OUTLINE 


• Computational domain is split into two regimes: 

□Non-linear source generation regime— Large-Scale Simulation. 

□Linear acoustic wave propagation regime— Surface-Integral Formulation (SIF). 

Linearized Euler. 

Kirchhoff ’s Method. 

Lighthill’s Theory. 
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Surface-Integral Formulation (SIF) 

Consider a cylinder of radius a, and length L enclosing the jet sound sources. 



x 



L 


H 


The mean flow outside this cylinder is stagnant and the disturbances are taken to be 
purely acoustic in nature, described by the simple wave equation 


c 8t 

The Fourier transform with respect to t is defined as 


p(m) 


| P(0 



Substituting in the wave equation, we obtain 

V^ + k^p = 0 k = — = 2nSM. S = 

c J Uj 

The integral solution is given by 

p(x o ,o>) = -j[G^-pf]dS 

where p(X Q , to) is the acoustic pressure at the observation point X 0 =(x 0 ^' 0 ,^ 0 ). G is 
the Green function, n is the noimal to the surface 5 and p(X, to) is the pressure dis- 
tribution on the surface at the point X=(x,r,<(>). 


NAS A/CP— 200 1-211152 923 


and 


3r J -» H m (qa) L 


(qa) q [ J 'm < qr)H m (qa > ~ J m (qa) H'm < 


For r=a, we can make use of the following relation 

2i 


H (u)J(u) -J (u)H(u) = 


TtU 


where u=qa. Substituting, we obtain 

m = ~ ~ M X o-jH ffi ( q r o ) 


^ = “T I e m COS [ m (*-*o)] J 

4rc am = 0 _« 

With G = 0 at the surface, the integral solution reduces to 

P = |p^adxd4> 

Substituting, we obtain the acoustic field as 

1 m = oo 

p ( V V *oj = — 2 /Jp (K ’ a - ♦> I e m cos [ m ( ♦ - ♦oil • 

4k m = 0 


J-oo H m (qa) x 


This is the final formula describing the relation between the acoustic far field and the 
pressure distribution on a cylindrical surface surrounding the jet noise sources* 
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Effect of grid spacing dr on pressure amplitude, L/a=30, dx=0.5. 
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OTHER TECHNIQUES FOR SOUND PROPGATION 

• Lighthilt Sound source may be noncompact 

• Kirchoof solution can be used to calculate the Fairfield sound given 
the pressure distribution over a cylindrical surface enclosing the 
source 

• An Analytical solution in terms of HankeL function can be used for 
extending the near field solution to the Fairfield- no pressure 
derivative term 
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THE SOURCE REGIME- VLES 


IMPROVED NUMERICAL SCHEME 
USE A COURSE MESH TO RESOLVE ONLY THE VLS 
MODEL THE UNRESOVED SCALES USING A HIGH-ORDER 
TURBULRENCE MODEL (K-epsilon) 
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Ve ry L a rge Eddy $ i m ulati on Me thod for Fast Prediction of Jet Ex ha u St Noise 


Centerline Axial Velocity Up Una Axial Velocity 



A* M Vttocfy Prom* at k/B - 229 



Axial Vetocty Promt afx/R ■ 6.32 



AH* Vaiocfty Promt *k/R- 10.45 
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r \ 

Vorticity Disturbance Magnitude 
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Instantaneous Pressure 
(VLES) 
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Technical Highlight 


3-D, Large-Eddy Simulation of a M=1.4 supersonic round jet excited by 
multiple frequencies at nozzle exit 



Mach Number 


Density 



Enhanced mixing and the breakdown of 
organized structures downstream of jet 
are observed. This near-field solution 
will be used to obtain far-field noise 
through the use of KirchhofF— type methods. 


Vorticity 
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Technical Highlight 

3-D, Large-Eddy Simulation of a M 
multiple frequencies at nozzle exit 



Vorticity 


1.4 supersonic round jet excited by 



Enhanced mixing and the breakdown of 
organized structures downstream of jet 
are observed. This near-field solution 
will be used to obtain far-field noise 
through the use of Kirchhoff-type methods. 
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NUMERICAL ISSUES 

• Extensive Computational Domain 60D x 40D 

• Helical modes- 3D 

• Various length scales: 

• Acoustic: Acoustic Wavelength 

• large-scale: diameter of the jet 

• Shear layer scales: momentum thickness 

• Mean flow 

• High-order scheme is needed to cover these scales in this large 
domain 

• Amplitude and phase accuracy (Dispersive and dissipation* errors 
need to be minimized) 

• Several Schemes are available for linear problems-Linear sound 
propagation can be handled. 

• Few Schemes are available for mixing noise (nonlinear): DRP, High- 
Order McCormick. Compact, Split Compact. 

• A suitable-scheme for acoustics-shock wave interaction is still being 
sought 

• Boundary Conditions still represent the major difficulty 
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OPEN PODIUM PRESENTATIONS 


Dr. J onathan Freund of UCLA presented results from DNS calculations. Flow field 
properties and far field noise were predicted. Comparisons made with the data of 
Stromberg for a low Reynolds number flow showed good agreement. 

Dr. Ray Hixon of NASA Glenn presented results from the application of a high 
order CAA code which used body-fitted curvilinear grids to compute nonlinear flows 
about complex geometries. Results for test cases from the 3rd CAA workshop were 
presented. Good results were obtained. 

Dr. Anastasios Lyrintzis from Purdue University presented an overview on Integral 
Methods for J et Aeroacoustics. Methods that can be used for this type of analysis 
include: Lighthi M's acoustic analogy, which requires a volume integral of the source 
terms; Kirchhoff's method, which requires a surface integral with a good flow solver 
to get a good solution on that surface; and the porous Ffowcs Williams - Hawkings 
equation. Dr. Lyrintzis recommended that work continue to improve refraction 
corrections and high frequency predictions. 

Dr. K. Viswanathan of Boeing presented an assessment of existing jet noise 
prediction methods. In particular, the MGBK method and the method of Tam were 
evaluated. A total of seven test cases, covering a range of jet Mach numbers and 
temperatures were used to evaluate these methods. Predicted spectra were 
compared with measured spectra for selected angles for each of the test cases. The 
same CFD solution (based on work of Thies and Tam) for the mean flow and 
turbulent energy was used for both noise prediction procedures. The method of Tam 
appeared to predict overall levels better than the MGBK method. Also, the MGBK 
method tended to under-predict the high frequency portion of the spectra while 
Tam's method predicted the spectral shape very well at angles away from the jet 
axis. Tam asserts that large scales are responsible for the noise near the axis and 
thus would not expect his method to predict the spectra shape at those angles. The 
MGBK method attempts to predict the spectra at all angles. In discussions after the 
presentation it was suggested that the better agreement of Tam's method with 
overall levels may be due to the fact that this method has been "calibrated” using 
aerodynamic input generated using the same method used in this study, whereas 
the MGBK code was "calibrated" using aerodynamic data from a different code. 
Though the choice of the CFD solution would play a role, as pointed out by 
Dr. Morris, there could be more fundamental issues with the acoustic analogy 
approach that would cause the observed discrepancy. 

Dr. Philip Morris of Penn State University presented information comparing the 
source terms in the Tam and MGBK methods. Dr. Morris' main conclusion was that 
the difference in predicted spectral shape between the two methods may be due to 
the frame of reference used to model the source auto-correlation function and the 
choice of model for the turbulent statistics. 


NAS A/CP— 200 1 -2 1 1 1 52 


939 



